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11 Conclusions 36 

Abstract 

In the present paper the linearized problem of plasma oscillations in layer (par- 
ticularly, in thin films) in external longitudinal alternating electric field is solved 
analytically. Specular - accommodative boundary conditions of electron reflection 
from the plasma boundary are considered. Coefficients of continuous and discrete 
spectra of the problem are found, and electron distribution function on the plasma 
boundary and electric field are expressed in explicit form. Absorption of energy of 
electric field in layer is calculated. 

Refs. 34. Figs. 2. 

Keywords: degenerate plasma, layer, specular-accommodative boundary condi- 
tion, plasma mode, expansion by eigen functions, singular integral equation, ab- 
sorption energy of electric field. 

PACS numbers: 52.35.-g, 52.90.+Z 



1 Introduction 

The present paper is devoted to degenerate electron plasma behaviour research. Anal- 
ysis of processes taking place in plasma under effect of external electric field, plasma waves 
oscillations with various types of conditions of electron reflection from the boundary has 
important significance today in connection with problems of such intensively developing 
fields as microelectronics and nanotechnologies [1] - [6]. 

The concept of " plasma" appeared in the works of Tonks and Langmuir for the first 
time (see [7]-[9]), the concept of "plasma frequency" was introduced in the same works 
and first questions of plasma oscillations were considered there. However, in these works 
equation for the electric field was considered separately from the kinetic equation. 

A. A. Vlasov [10] for the first time introduced the concept of "self-consistent electric 
field" and added the corresponding item to the kinetic equation. Now equations describing 
plasma behaviour consist of anchor system of equations of Maxwell and Boltzmann. The 
problem of electron plasma oscillations was considered by A. A. Vlasov [10] by means of 
solution of the kinetic equation which included self- consistent electric field. 

L.D. Landau [11] had supposed that outside of the half-space containing degenerate 
plasma external electromagnetic field causing oscillations in plasma is situated. By this 
Landau has formulated a boundary condition on the plasma boundary. After that the 
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problem of plasma oscillation turns out to be formulated correctly as a boundary value 
problem of mathematical physics. 

In [11] L.D. Landau has solved analytically by Fourier series the problem of coUisionless 
plasma behaviour in a half-space, situated in external longitudinal (perpendicular to the 
surface) electric field, in conditions of specular reflection of electrons from the boundary. 

Further the problem of electron plasma oscillations was considered by many authors. 
Full analytical solution of the problem is given in the works [12] and [13] . 

This problem has important significance in the theory of plasma (see, for instance, [2], 
[14] and the references in these works, and also [15], [16]). 

The problem of plasma oscillations with diffuse boundary condition was considered in 
the works [17], [18] by method of integral transformations. In the works [19], [20] general 
asymptotic analysis of electric field behaviour at the large distance from the surface was 
carried out. In the work [19] particular significance of plasma behaviour analysis close 
to plasma resonance was shown. And in the same work [19] it was stated that plasma 
behaviour in this case for conditions of specular and diffuse electron scattering on the 
surface differs substantially. 

In the works [22] and [22] general questions of this problem solvability were considered, 
but diffuse boundary conditions were taken into account. In the work [22] structure of 
discrete spectrum in dependence of parameters of the problem was analyzed. The detailed 
analysis of the solution in general case in the works mentioned above hasn't been carried 
out considering the complex character of this solution. 

The present work is a continuation of electron plasma behaviour in external longitu- 
dinal alternating electric field research [22] - [26] . 

In the present paper the linearized problem of plasma oscillations in external alter- 
nating electric field in layer (particularly, in thin films) is solved analytically. Specular 
- accommodative boundary conditions for electron refiection from the boundary are con- 
sidered. In [24]-[26] diffuse boundary conditions were considered. 

The coefficients of continuous and discrete spectra of the problem are obtained in 
the present work, which allows us to derive expressions for electron distribution function 
at the boundary of conductive medium and electric field in explicit form, to reveal the 
dependence of this expressions on normal momentum accommodation coefficient and to 
show that in the case when normal electron momentum accommodation coefficient equals 
to zero electron distribution function and electric field are expressed by known formulas 
obtained earlier in [12], [13]. 

The present work is a continuation of our work [27] , in which questions of plasma waves 
specular refiection from the plane boundary bounding degenerate plasma were considered. 

Let us note, that questions of plasma oscillations are also considered in nonlinear 
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statement (see, for instance, the work [28], [29]). 



2 Formulation of problem 

Let degenerate plasma occupy a slab (particularly, thin films) —a < x < a. 
We take system of equations describing plasma behaviour. As a kinetic equation we 
take Boltzmann — Vlasov r-model kinetic equation: 

^ + v|^ + eE|i = M£dW(£lM. (1.1) 
at or dp T 

Here / = /(r, v,t) is the electron distribution function, e is the electron charge, 
p = mv is the momentum of an electron, m is the electron mass, r is the character time 
between two collisions, E = E(r, t) is the self-consistent electric field inside plasma, feq = 
feqir, t) is the local equilibrium Fermi — Dirac distribution function, f^q = Q{EF{t, x) — ^), 
where Q{x) is the function of Heaviside, 



1, X > 0, 
0, X <0, 



8.p(t,x) = ^mvp(t,x) is the disturbed kinetic energy of Fermi, £ = ^mv^ is the kinetic 
energy of electron. 

Let us take the Maxwell equation for electric field 

divE(r,i) = 47rp(r,t). (1.2) 

Here p(r, t) is the charge density, 

p(r, t) = ej (/(r, v, t) - /o(v)) dO^, (1.3) 

where 

d^F = ^^^^3 , d^P = dp^dpydp^. 
Here /o is the undisturbed Fermi — Dirac electron distribution function, 

/o(£) = e(£^-£), 

h is the Planck's constant, u is the effective frequency of electron coUisions, u — l/r, 
Ep — ^mvp is the undisturbed kinetic energy of Fermi, vp is the electron velocity at the 
Fermi surface, which is supposed to be spherical. 
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We assume that external electric field outside the plasma is perpendicular to the 
plasma boundary and changes according to the following law: Eq exp{—iujt) . 

Then one can consider that self-consistent electric field E(r, inside plasma has one 
x-component and changes only lengthwise the axis x: 

E^{E,{x,t),0,0}. 

Under this configuration the electric field is perpendicular to the boundary of plasma, 
which is situated in the plane x — 0. 

We will linearize the local equilibrium Fermi — Dirac distribution fgg in regard to the 
undisturbed distribution /o(£): 

/e, = /o(£) + [£F(a;,t)-£]5(£F-£), 

where S{x) is the delta - function of Dirac. 

We also linearize the electron distribution function / in terms of absolute Fermi — 
Dirac distribution /o(£): 

/ = /o(£) + /i(^,v,i). (1.4) 

After the linearization of the equations (1.1)-(1.3) with the help of (1.4) we obtain 
the following system of equations: 

^ + + uf,(x, V, t) = 5(8,F - £) [eE,(x, t)v, + v[Ef{x, t) - Ep]] , (1.5) 



dt dx 

dEx{x,t) 87re 



j Mxy,t)d'p' (1.6) 



we find: 



dx {2TiKf 
From the law of preservation of number of particles 

[Zf(x, t) - £f] j SiEp - e)d'p = j hd?p. (1.7) 

From the equation (1.5) it is seen that we should search for the function /i in the form 
proportional to the delta - function: 

fi^5(8,F-8,)H(x,ii,t), = (1.8) 

V 

The system of equations (1.5) and (1.6) with the help of (1.7) and (1.8) can be trans- 
formed to the following form: 

OH OH V 

'of ^ ^^^~0x ^ ^ evFnE^{x, t) + 2 J ^' ' 
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— ' H[x,n,t)an. 



dx {2nhy 7_i 

Further we introduce dimensionless functions 

E,{x,t) H{,v.jLt) 

and pass to dimensionless coordinate xi — x/a. We obtain the following system of 
equations 

dH dH 1 

+ H-^ + i^H{xi,iJ,,ti) ^ He{xi,ti) + - J H{xi,n',ti)dn', (1.9) 

de{xi,ti) IQ-n'^e'^m'^VFa^ '"^ 



j H{xi,fi',ti)dfi'. (1.10) 



9X1 (27r;i)3 

Here cUp is the electron (Langmuir) frequency of plasma oscillations, 

^ m 

A?^ is the numerical density (concentration), m is the electron mass. 

We used the following well-known relation for degenerate plasma for the conclusion of 
the equations (1.9) and (1.10) 

h / 



3 Boundary conditions statement 

Let us outline the time variable of the functions H{xi, fj,,ti) and e{xi,ti), assuming 

H{xun,h)^e-"'^'^h{xuii), (2.1) 

e(a:i,ti) = e-*'^i*ie(a:i). (2.2) 

The system of equations (1.9) and (1.10) in this case will be transformed to the fol- 
lowing form: 

1 

dfx f 

H-Q^ + WQh{xi, jj) = jjLe{xi) + j h{xi, jji')djji\ (2.3) 



-1 

1 



dxi 2 



h{xi,n)dn. (2.4) 
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Here 

d , . s av auj 2 nfd'^v\'^ 

Vf Vf Vf \ vf ^ 

The constant Up can be expressed through Debaye radius td 

2 Vf 
u = 3-5-, td = — ■ 

Further instead of Xi we write x. We rewrite the system of equations (2.3) and (2.4) 
in the form: 

1 

IJ,^ + woh{x,iJ,) = iJ,e{x) + ^ J h{x,ii)dn' . (2.5) 



de{x) 



^ y h{x,i,')di,'. (2.6) 



dx 

For electric field in plasma on its border the boundary condition is satisfied 

e(-l) = e„ e(+l) = e„ (2.7) 

where 1 is the dimensionless depth (width) of a semilayer. 

Condition of symmetry of boundary conditions (2.7) and the equations (2.5) and (2.6) 
mean, that electric field e{x) in the layer possess properties of symmetry 

e{x) = e{-x). (2.7a) 

The non-flowing condition for the particle (electric current) flow through the plasma 
boundary means that 

1 1 
J lih{-l,ii)dii^ J lih{l,ii)dii^O. (2.8) 
-1 -1 

In the kinetic theory for the description of the surface properties the accommodation 
coefficients are used often. Tangential momentum and energy accommodation coefficients 
are the most-used. For the problem considered the normal electron momentum accom- 
modation under the scattering on the surface has the most important significance. 

Owing to properties of symmetry of electric field and distribution function concerning 
a plane — middle of a layer — further let's enter the coefficient of accommodation of a 
normal momentum through momentim of electron streams on the bottom surface of the 
layer. 
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The normal momentum accommodation coefficient is defined by the foUowing relation 

«p = O^c^p^h (2.9) 

where Pi and Pr are are the flows of normal to the surface momentum of incoming to the 
boundary and reflected from it electrons, 



Pi^ J n^h{-l,n)dfx, (2.10) 



-1 



Pr^ J fj,'^h{-l,fj,)dfx, (2.11) 







quantity Pg is the normal momentum flow for electrons reflected from the surface which 
are in thermodynamic equilibrium with the wall, 



Ps 



1 

J iJ.X(ii)di^, (2.12) 



where the function 

hsin) =As, < < 1, 

is the equilibrium distribution function of the corresponding electrons. This function is 
to satisfy the condition similar to the non-flowing condition 

1 

/x/i(— 1, /x)ci// + / iJihs{ij)diJL — Q. (2-13) 



-1 



We are going to consider the relation between the normal momentum accommodation 
coefficient and the diffuseness coefficient q for the case of specular and diffuse boundary 
conditions which are written in the following form 

/i(-l,/x) = (l-g)/i(-l,-/x) + a„ 0</x<l. (2.14) 

Here q is the diffusivity coefficient (0 ^ ^ 1), is the quantity determined from the 
non-flowing condition. 

Prom the non-flowing condition we derive 

10 1 

j lih{—l,ii)dii = J iJ,h{—l, iJ,)dfi + J iJ,h{—l, iJ,)diJ, ~ 0. 

-1 -1 
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In the second integral we replace the integrand according to the right-hand side of 
the specular-diffuse boundary condition (2.14). After that, using the obvious change of 
integration variable, we obtain that 



as — —'2q J fj,h{—l, fj,)dfj,. 
-1 

Let us use the boundary condition (2.13). Using the analogous to the preceded hne of 
reasoning we get 



Ag — —2 J fj,h{—l, fj,)dfj,. 



From the two last equations we find that 

as = qAs. (2.15) 
Further we find the difference between two flows 

1 

Pi — Pr = J iJ?h{—l^ jji)djji — J jj?h{—l^ jji)diJi. 

-1 

In the second integral we use the boundary condition (2.14) again. With the help of 
(2.15) we obtain that 

1 



Pi — Pr — q J fj,^h{—l, fj,)dfj, — J fx^agdfj, — 



^q J /x^/i(-l, l^jdji -q J ji'^Asdji = qPi - qPg. 

-1 

Substituting the expressions obtained to the definition of the normal momentum ac- 
commodation coefficient, we have 

- Pr qPi - qPs 
ttp = = = q. 

Pi Ps Pi Ps 

Thus, for spccTilar - diffuse boundary conditions normal momentum accommodation 
coefficient ap coincides with the diffuseness coefficient q. 

Equally with the specular - diffuse boundary conditions another variants of boundary 
conditions are used in kinetic theory as well. 

In particular, accommodation boundary conditions are used widely. They are divided 
into two types: diffuse - accommodative and specular - accommodative boundary condi- 
tions (see [30]). 
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We consider specular - accommodative boundary conditions. For the function h this 
conditions will be written in the following form 

= + Ao + Ai/x, 0</x<l. (2.16) 

/i(l,//) = /i(l,-//) - Ao + Ai//, -l<//<0. (2.17) 

Let's notice, that (2.17) the same appearance, as (2.16) has conditions. Really, let's 
replace on — in the condition (2.17). We will receive a condition 

h{l,-jji) = h{l,jjL) - Aq- Aijji, < < 1, 

whence we receive in accuracy a condition (2.16). It means, that both conditions (2.16) 
and (2.17) it is possible to write down in the form of one 

h{±l,ii) = h{±l,-ii) + Aq + Ain, < // < 1. 

If in (2.16) we assume Aq = Ai = 0, then specular - accommodative boundary condi- 
tions pass into pure specular boundary conditions. 

Coefficients and Ai can be derived from the non-flowing condition and the definition 
of the normal electron momentum accommodation coefficient. 

The problem statement is completed. Now the problem consists in finding of such 
solution of the system of equations (2.5) and (2.6), which satisfies the boundary conditions 
(2.7) and (2.16). Further, with the use of the solution of the problem, it is required to built 
the profiles of the distribution function of the electrons moving to the plasma surface, and 
profile of the electric field. 



4 The relation between flows and boundary 

conditions 

First of all let us find expression which relates the constants Aq, Ai from the boundary 
condition (2.16). To carry this out we will use the condition of non- flowing (2.12) of the 
particle flow through the plasma boundary, which we will write as a sum of two flows 

1 

Nq = J fj,h{—l, fj,)dfj, + J fj,h{—l,fj,)dfj, — 0. 

-1 

After evident substitution of the variable in the second integral we obtain 

1 





Io = J II /i(-l,//) - /i(-l, -//)jd// = 0. 
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Taking into account the relation (2.16), we obtain that = 

With the help of this relation we can rewrite the condition (2.16) in the following form 

h{-l, ii) - M-l, -/.) = Ai(/. - ^), < /X < 1, (3.1) 

or, that is equivalent, 

/i(l, ii) - /i(l, -/.) = yli(/. - ^), < < 1, (3.1') 

We consider the momentum flow of the electrons which are moving to the boundary. 
According to (3.1) we have 

Pi^Pr- (3.2) 

It is easy to see further that 

Ps = y . (3.3) 

With the help of the formulas (3.2) and (3.3) we will rewrite the definition of the 
accommodation coefficient (2.9) in the form 

As Ai 

apPr- + —{I- ap)^Q. (3.4) 
Let us consider the condition (2.13). We rewrite it in the following form 



— ^ + / /i/i(— 1, /i.)(i/x = 0. 



-1 

From this condition we obtain 





Ag — —2 J fj,h{—l, iJ,)diJ, = 2 J nh{—l,—iJ,)diJ,. 

-1 

Using the condition (3.1), we then get 

1 

As ^2 iih{-l,ii)dii. (3.5) 



Now with the help of the second equality from (2.11) and (3.5) we rewrite the relation 
(3.4) in the integral form 

1 

«P y (/^^ - ^/^)^(-l>/^)c^/^ = -^(1 - «p)A- (3.6) 



Now the boundary problem consists of the equations (2.5) and (2.6) and boundary 
conditions (2.7), (3.1) and (3.6). 
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5 Separation of variables and characteristic system 

Application of the general Fourier method of the separation of variables in several 
steps results in the following substitution [31] 

hr,{x, ii) = exp(-^)$(T7, /x) + exp(— )<l>(-r?, /x), (4.1) 

T) T) 

e,{x) = [exp(-^) +exp(^)]E(r/), (4.2) 

where rj is the spectrum parameter or the parameter of separation, which is complex in 
general. 

We substitute the equalities (4.1) and (4.2) into the equations (2.5) and (2.6). We 
obtain the following characteristic system of equations 

1 

wo{ri - /x)$(77, //) = wE{v) + \j *(?7, l^W. (4.3) 



-1 
1 



w^{r] + /x)$(-?7, /i) = wEiri) + <!>(-??, fi')dfj,', (4.4) 



-1 

1 



-1 
1 



-1 

Prom the equations (4.5) and (4.6) we obtain 

1 1 
j $(77, n)dn = - y" $(-77, n)dii. (4.7) 
-1 -1 

Let us introduce the designations 

1 

7i(77) = j ^v,l^)dl^- (4.8) 
-1 

Prom the equation (4.5) we find, that 

E{v) = -TT^Vniv), (4.9) 



2w, 
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whence 

Wo E{r)) 

n[V) = -2— . 

u; V 

By means of equalities (4.7) - (4.9) we will copy the equations (4.3) and (4.4) 

{rj - /x)$(r;, II) = ^{fiv - v!): (4.10) 

Wo 

{rj + /x)$(-r], /x) = ^{i,^ + nl). (4.11) 

Wo 

Here 

2 t'o'ii'o i^(i^ — iu;) £^^0 ^ ^ .oj ^ .VL 

?7i = = — — = ^ = — ' 2;o = 1 -^- = 1 12 = —. 

Solution of the system (4.10) and (4.11) depends essentially on the condition whether 
the spectrum parameter 77 belongs to the interval — 1 < 77 < 1. In connection with this 
the interval — 1 < 77 < 1 we will call as continuous spectrum of the characteristic system. 

Let the parameter 77 e (—1,1). Then from the equations (4.10) and (4.11) in the class 
of general functions we will find eigenfunction corresponding to the continuous spectrum 

$(77, /.) = EMp^^l^ + g^^ri)5{ri - Ai), (4.12) 
zo V ~ 1^ 

$(-77, /.) = ^pEl±li + g,(rj)S{rj - /.). (4.13) 

Zo V + 

In these equations (4.12) and (4.13) 6{x) is the delta-function of Dirac, the symbol 
Px~^ means the principal value of the integral under integrating of the expression x~^. 

Substituting now (4.12) and (4.13) in the equations (4.5) and (4.6), we receive the 
equations from which we obtain 

= -2^^£(,). = = 2^M£(,). (4.14) 

7/ '^p V 

Here dispersion function is entered 

1 

A(.) = 1 - f / ^^^^dn, (4.15) 
2c J 11 — z 
-1 

where ^ ^ ^ 

Wq e Zq 2 



ul 3 



Functions (4.12) and (4.13) are called eigen functions of the continuous spectrum, 
since the spectrum parameter 77 fills out the continuum (— 1,+1) compactly. The eigen 
solutions of the given problem can be found from the equalities (4.1) and (4.2). 
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Substituting relations (4.14) in (4.12) and (4.13), we will present last expressions in 
the following form 



$(r^,//) = 
$(-77,//) - 



Wo 

_ Ejy) 

Wo 



ri- n 



V -I 



or 



$(77, /.) = ^F(77, /x), $(-77, /.) = ^F(-7^, /x). 



Wo 



Wo 



where 



7;-/x t; 



Xv) 



r] + n T) 



It will be necessary for us the following relation of symmetry F{—r], —fx) = —F[r], jji). 
Let us notice, that eigen function F{ri, n) satisfies to following condition of normaliza- 



tion 



/ F{±ri,fj,)dn = T2cP-. 
J V 



So, eigen function of a continuous spectrum is constructed and it is defined by equality 



XWo 



hr^{x,n) = [exp - ^jF(7/,/x) + exp (^^^F{-ri,n) 
or, in explicit form. 



Wo 



XWo\ HT] - Tjf fXWo\fiV + Vf 

exp h exp 

V n / n — a \ n / 



r) / rj + jj, 



-2c 



exp ^^5(77-//) - exp (^^) 5(77 + /i) 



Wo 



Let us replace exponents by hyperbolic functions 

fXWo\ ,XWo , ,XWo ( XWo\ , 

exp ( ) = ch h sh , exp I I = ch 

V77/ 77 7? \ n / 



XWo 1 XWo 

sh 



77/ 77 77 \ r] / 77 77 

and also we will transform both square brackets from the previous expression. As a result 
we receive, that for the first square bracket it is had 



exp 



XWo \ 117] — 7]l 
T] ) T] — jJi 



+ exp 



/ xwo \ fir] + r]( 



15 



= exp 



xwo\ (iiri - 7]f){7] + n) ^ ^^^^ /xwo\{ix7] + r]f){r]-ii) 
J rf — (J? \ T) J 



V 



77^ — n'^ 



= 2 



rj A* ?7 
For the second square bracket we have 



77 ^"^ — rf 



= - 5{ri- 11) + 5{ri + n) 
As a result we receive, that 



exp - exp (p^^5{ri + n) = 



sh h 



ch 



XWq 



hJx,ij) — ^< eh ^ — - 

Wo y 7] 



\ T] - n V + 1^ 



+2c^{-5{r, - f,) + 5{rj + I,)) 



+ sh 



XWq 

V 



+2c^{S{ri-ii) + S{r, + ii)) 



We will designate further 

\ n — u n + u / — 



and 



^{r],ii)^-F{r],ii) + F{-r],ii)^P 



+ — ) = 2P- 

II rj + n J 7]^ — 11^ 



Thus, eigen function of a continuous spectrum it is possible to present in the form of 
a linear combination of a hyperbolic sine and kosine 



hjj{x,n) — ch 



Wo 



rf — iP 



- c {^{ri -11)- S{ri + 11)) 



V 



sh 



XWo 

V 



rf — iJ? 



V 



Let us notice, that eigen functions of a continuous spectrum it is possible to present 
and in such form 



hr^{x,lj) = 



Wo 



ch ^ [Fin, ii)+F{-n, 11)) +sh ^ [-Fin, ii)+F{-n, 11) 



. (4.16) 
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The dispersion function X{z) we express in the terms of the Case dispersion function 

[31] 



where 



1 1 

1 f T dr 



T — Z 



-1 



is the Case dispersion function. 

In the complex plane dispersion Case function is calculated through the logarithm 



Ac(^) = l + ^ln^, zeC\[-l,l], 
z z + 1 



and on cut by formula 



Ao(r7) = l + |lni-p^, 77 e (-1,1). 



The boundary values of the dispersion function from above and below the cut (interval 
(— 1, 1)) we define in the following way 

A^(/x)= lim A(/x±z£), /xe(-l,l). 

The boundary values of the dispersion function from above and below the cut are 
calculated according to the Sokhotzky formulas 



from where 



X^(^i)^X(^i)±-^(rjt-^i% -1</.<1, 

ZTjiZo 



A+(/.)-A-(/.) = ^M^?-/^'), ^^^^^^^^ = A(/.), -1</.<1, 

where 

1 

A(/.) = l + -^ / ^^^dr), 

-1 

and the integral in this equality is understood as singular in terms of the principal value 
by Cauchy. Besides that, the function A(/x) can be represented in the following form 

A(/.) = 1 - 1 + 1 (l - ^) Ao(/.), Ao(/.) = 1 + f In 1^. 
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6 Eigen functions of discrete spectrum and plasma 

waves 

According to the definition, the discrete spectrum of the characteristic equation is a 
set of zeroes of the dispersion equation 

M.o. (6.1) 

z 

We start to search zeroes of the equation (5.1). Let us take Laurent series of the 
dispersion function 

A(^) = Aoo + ^ + ^ + ---, \z\>\. (5.2) 



Z" ^4 



Here 



^ _ ^ , , 1 1 ^ 1 /I 1 \ ^ 1/11 

A.. = A(oo) = l-- + ^, A. = --(--^j, A4 = --(--^ 

We express these parameters through the parameters 7 and e 

2(Q - 1) + i£ + (Q - 1)(Q - 1 + is) 



Aoo = A(oo) 



(Q + iz) 



2 



_ ^^hie{i\^ie) 15 + 7i£:(fi + «£:) 

Ao — _ — t^tt: — , A4 



15(Q + i£)2 ' " 35(Q + i£)2 ■ 

It is easy seen that the dispersion function (4.9) in coUisional plasma (i.e. when £ > 0) 
in the infinity has the value which doesn't equal to zero: Aoo = A(oo) 7^ 0. 

Hence, the dispersion equation has infinity as a zero 77^ = 00, to which the discrete 
eigensolutions of the given system correspond 

Wo 

This solution is naturally called as mode of Drudc. It describes the volume conduc- 
tivity of metal, considered by Drude (see, for example, [32]). 

Let us consider the question of the plasma mode existence in details. We find finite 
complex zeroes of the dispersion function. We use the principle of argument. We take 
the contour (see Fig. 1) r+ = Fr U 7^, Tr = {z : \z\ = R, R — 1/e, e > 0}, which 
is passed in the positive direction the cut [—1,-1-1], and which bounds the biconnected 
domain Dr. 

Let us notice, that dispersion function in area Dr has no poles. Then owing to the 
principle of argument [33] zeroes number N in area Dr it is equal 




Fig 1. Contour = Fij U for calculations of number of zero 
dispersion function. 
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Considering the limit in this equahty when £ — )■ and taking into account that the 
dispersion function is analytic in the neighbourhood of the infinity, we obtain that 

111 



N= — I d InA+fr) -—Id In A-(t) = — I din 
-1 -1 -1 

So, we have received, that 



1 A+(r) 



= — / ci In 



2niJ A-(t)' 

-1 

We divide this integral into two integrals by segments [—1,0] and [0, 1]. In the first 
integral by the segment [—1,0] we carry out replacement of variable r —r. Taking into 
account that A"'"(— r) = A~(r), we obtain that 



1 1 



N ^ — din ; ; = — din _ { = - arg — ^ ^ 



27ri J X (r) ni J A (r) tt A (r) 



1 



(5.3) 



A+(r) 

Here under symbol arg G(r) = arg we understand the regular branch of the 

A (r) 

argument, fixed in zero with the condition: arg (7(0) = 0. 

We consider the curve = Ag'(7,£) : z = G{t), ^ t ^ +1, where 

A+(r) 



G(t) 



X-(t) 



It is obvious that G(0) = 1, lim G(t) — 1. Consequently, according to (5.3), the 

T—^ + l 

number of values A'^ equals to doubled number of turns of the curve 7 around the point 
of origin, i.e. 

N^2h{G), (5.3') 

where x(G) = Ind[o,+i]G(T) is the index of the function G{t). 

Thus, the number of zeroes of the dispersion function, which are situated in complex 
plane outside of the segment [—1,1] of the real axis, equals to doubled index of the function 
G(t), calculated on the "semi-segment" [0, +1]. 

Let us single real and imaginary parts of the function G(/i) out. At first, we represent 
the function (j(/x) in the form 

^/ X _ (^0 - l)r]l + {rji - //^)Ao(/x) + is{ii){r]l - /x^) 

"^^^^ (^0 - iH + {vi - /^^)Ao(/^) - isiiJ^m - /^^) ■ 

where ^ 

= A(/x) = 1 h — (l-^j Ao(/i), 



20 



and 



Ao(/.) = l + ^ln.^ ^ 



2 l + 

is the dispersion function of Case, calculated on the cut (i.e., in the interval (—1, 1)). 
Taking into account that 



we obtain 



where 



V £' " 3 3 3' vu,,, 3 3 

p-(/i)+^Q-(/i) 



p+(/x) + ^g+(/i) 



g±(Ai) = £0(1 + Ao(Ai)) ± s{^J){s'' - 3fi^). 
Now we can easily separate real and imaginary parts of the function G{ii) 

Here 

^(/x) = [p+(^)]2 + [g+(/.)]2 = 

+[n^ + Ao(3/i2 - - £Qs]^ + [en(l + Ao) - 5(3//^ - e^)]^, 
giifi) = P^{fx)P-{fi) + Q+(/x)g-(/i) = 
= [Q' + Ao(3/.^ - - £202[s^ - (1 + Ao)^] - (3/.^ - e'fs', 

+2s[Q2(3//2 - + Ao(3//2 - £2)2 + £2q2(i + Ao)], 

We consider (see Fig. 2) the curve L, which is defined in implicit form by the following 
parametric equations 

L = {(Q,£): g'i(//;Q,£) = 0, g'2(Af;f^,£) = 0, ^ // ^ l}, 

and which lays in the plane of the parameters of the problem (7,£), and when passing 
through this curve the index of the function G(/i) at the positive "semi-segment" [0, 1] 
changes stepwise. 

From the equarion = we find 

- 3/.2 + £2Ao(/.) ■ ^'-^^ 
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Now from the equation gi — with the help of (5.4) we find that 

e = \fL^, 

where 

L2{ix) = 



(5.5) 



3;u2s2(/i) 



Ao(/x)[.2(;.) + (l + Ao(/.))2] 
Substituting (5.5) into (5.4), we obtain 



where 



Q = +VLi(At), 

3/.2[52(;,) + Ao(/.)(l + Ao(/.))]^ 



(5.6) 



Ao(//)[s2(/.) + (l + Ao(/i))2] ■ 
Functions (5.5) and (5.6) determine the curve L which is the border if the domain 
(we designate the external area to the domain as D~) in explicit parametrical form 
(see Fig. 2). As in the work [34] we can prove that if (7, s) G D^, then x(G) = 
Ind[o,+i]G(/x) = 1 (the curve L encircles the point of origin once), and if (7, e) e then 
x(G') = Ind[o,+i]G(//) = (the curve L doesn't encircle the point of origin). 

We note, that in the work [34] the method of analysis of boundary regime when 
(7, e) e L was developed. 

From the expression (3.2) one can see that the number of zeroes of the dispersion 
function equals to two if (7, e) e D+, and equals to zero if (7, e) € D~. 

Since the dispersion function is even its zeroes differ from each other by sign. We 
designate these zeroes as following ±770, by rjo we take the zero which satisfies the condition 
Re r]o > 0. 

The following solution corresponds to the zero ±770 



hr,o{x,ll) 



Here 



ch^[F(?7o,/i) 
^0 



F{-r]o,n)] +sh 
er,o{x) = 2ch 



XWo 

Vo 



[-F{r]o,fi) + F{-r]o,^)], 



(5.7) 
(5.8) 



W - Vt 



It is easy to see, that function /ir?o(^) A*) even on rjo 
Function h^^{x,^) we will present in the explicit form 



?7o/^ + m 

110 + fJ' 



hr,o{x,IJ') = ch 



Vo 



sh 



XWq 

Vo 



(5.9) 




Fig. 2. 
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or 



n) = ch (/)(?7o, li) + sh V(?7o, A*)- (5.9 ) 

^0 ?7o 



Here 



^^0, /^) = -i^ r^o, + F{-rio, IJ) = + — = 2 —■ 

This solution is naturally called as mode of Debay (this is plasma mode). In the 
case of low frequencies it describes well-known screening of Debay [3] . The external field 
penetrates into plasma on the depth of r^, is the raduis of Debay. When the external 
field frequencies are close to Langmuir frequencies, the mode of Debay describes plasma 
oscillations (see, for instance, [3, 32]). 

Note 5.1. If to enter expression c/z "inside" of expression of dispersion function we 
will receive expression for dispersion function h{z) from our article [21] 



7 Expansions by eigen functions 



We will seek for the solution of the system of equations (2.5) and (2.6) with boundary 
conditions (3.1), (3.6) and (2.7) in the form of linear combination of discrete eigen solutions 
of the characteristic system and integral taken over continuous spectrum of the system. 
Let us prove that the following theorem is true. 

Theorem 6.1. System of equations (2.5) and (2.6) with boundary conditions (3.1), 
(3.6) and (2.7) has a unique solution, which can he presented as an expansion by eigen 
functions of the characteristic system 



h{x,ii) 



Wo Wo 



— ^exp(- 
-/^ 



WoX 

Vo 



WqX 

Vo 



+ 



1 r 



+ 



-1 



exp(- ^)F(77,/x) + exp (^)F(-7y,/x) 



e{x) ^ + Eo exp ( ) + exp ( j 

L T] ^ Vo 



rWoX-, 



E{v) 

Wo 

+ 



dr], 



(6.1) 
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+ 



exp ( ) + exp ( ) 



,WqX, 



E{r)) drj. 



(6.2) 



Here Eq and is unknown coefficients corresponding to the discrete spectrum (Eq is 
the amplitude of Debay, Ei is the amplitude of Drude), E[rj) is unknown function, which 
is called as coefficient of continuous spectrum. 

When {fl,e) G D~ in expansions (6.1) and (6.2) we should take Eq = 0. 

Further we will consider the following case {^,e) G D~^. 

Our purpose is to find the coefiicient of the continuous spectrum, coefficients of the 
discrete spectrum and to built expressions for electron distribution function at the plasma 
surface and electric field. 

Proof. Let us notice, that the formula (6.1) can be transformed and to such form 



h{x,ii) 



Eoo , 2Eq 



Wo Wo 



r]o{fi^ - vf ^-^ wox _ iJ,{r]l - r]j) wpx 



+ 



+ 



^1 

Wo J 



-1 



eh '^^^F(ri, n) + sh ^^^^F(— 77, /x) 
77 77 



E{ri)dri. 



(6.10 



Let us consider expansion (6.1), we will replace in it // on — Then we will sub- 
stitute the difference h{l,fj,) — h{l, —/i) in a boundary condition (3.1). After variety of 
transformations let us have 



EoofJ' + Eo 



FivQ,l^) + F{-r)o,n) 

zoA 



ch^+ 
Vo 



+ y"F(r?,/x)E(r?)ch^rfr?=^(/x-^), < /x < 1. 



Substituting expansion (6.2) in (2.7), we will have 



(6.3) 



1 

£^00 + 2^0 ch — + 2 / E(ri) ch — dri = 1. 
Vo J V 



(6.4) 



Let us pass from Fredholm integral equation (6.3) to singular integral equation with 
Cauchy kernel, having substituted in (6.3) obvious representation F{rj,ii) 

1 



E^lji + Eoip(r]o, //) ch — + 

Vo 



- ^^.£;(77) ch ^rf7,-2cM eh ^ £;(;.) 



r]- ^ 



V 



WqA 



< n <1, 



(6.5) 
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where 

ip{r)Q, = F 770, ^J) + F[-r]o, n) = + ■ = 2/1^ -. 

Vo-I^ Vo+I^ Vo-I^ 

It is easy to check up, that function 

1 ^ 

M(z) ^ / ~ ^1 eh — EMdri (6.6) 
J r] — z V 

-1 

is odd. Besides, all members of the equation (6.5) are odd on /i, except one member from 

its right part —-wqAi. 

Hence, the equation (6.5) can be extended in the interval — 1 < // < 1 in the next 
symmetric form 

1 ^ 
E^fi + Eoip{r]o,n)ch. — + / ch — ^^^^ — —E{rf)drj- 

-1 

-2cM ch ^ Eit,) - = -\woA, sign -1< < 1. (6.7) 



Let us reduce the equation (6.7) to Riemann — Hilbert boundary value problem. For 
this purpose we will take advantage Sohkotsky formulas for the auxiliary functions M[z) 
and dispersion function A(^) 

M+{^i)- M-{ii)^2mch — {ii'^ -r]l)E{ii), -!<//<!, (6.8) 

A* 

where 

M(;.)= /^^ch^F(,)ci,, 
-1 

besides, last integral is understood as singular in the sense of principal value of Cauchy, 
and 

A+(/.)-A-(/.) = -^M/^^-ryD, ^^^^^^±^ = A(/.), -!</.<!, 
c z 

where 

A(,.) = 1 - f / ^d,'. 

2cJ n' - 1^ 



26 



besides, last integral is understood as singular in the sense of principal value of Cauchy 
also. 

As result of use of last formulas we will come to the boundary value problem 



woA 



1 1 . 
2 -At = --wqAi sign II, -!<//<!. 

We transform this equation to the form 

^(M+ + M-)(A+ - A") + ^(M+ - M-)(A+ + A-)+ 

+(A+ - A") (e^ii + Eo(p{r]o, /x) ch — - ^^^^lA = 

1 tlT 

= -WqAi sign - r]l), -1 < < 1. 

6 C 

From here we receive the following boundary condition of boundary value Riemann 
Hilbert 

A+(/x) \m+{ij,) + E^fx + Eoifiivo, /x) ch ^ - J - 
= A"(/x) [m-(/x) + E^ii + Eo(p{r)o, /x) ch — - ^^^^/j] = 



m 
3c 



-1 < < 1. 



(6.9) 



We rewrite this problem in the form 



- $ (/i) = ^woAi/x(/x^ - ?7i) sign n, -1 < < 1. 



(6.10) 



In the problem (6.10) $ (//) are boundary values on interval — 1 < < 1 of function 



^(z) = X+(z) 



M+{z) + E^z + Eo(^(77o, ^) ch ^ - 



]■ 



which is analytic in complex plane with cut C \ [—1,1]. 

The problem (6.10) is a problem of special case about jump. The jump problem is the 
problem of finding of analytical function by its jump on the contour L: 

$+(/.)-*+(/.) = (^(/.), pteL. 

The solution of such problems in a class of functions decreasing at infinitely remote 
point it is given by integral of Caucy type 

ip{T)dT 



^ ^ 2m J T 
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V?(/i) = ^WoAindj,^ - rjl) sign /i, -1 < < 1, 



However, ih the problem (6.10), when 

2m 
6c 

unknown function ^{z) has at infinitely remote point the following asymptotic 

^{z) = 0{z), z ^ oo. 

Therefore it is necessary to search the solution of the problem (6.11) in a class of the 
growing as z in the vicinity of infinitely remote point. 

According to [33] the general solution of problem (6.10) it is given by the formula 

ip{T)dT 



^ ' 2m J T 



+ Ciz. 



In explicit form the general solution of problem (6.10) write down the following form 
X(z) M(z) + Eoip(r]o, 2;) ch — + E^wq - 



where Ci is the arbitrary constant. 
Let us enter auxiliary function 



/i(/i 



//J j sign fj, 



fjL — Z 



T{z) 



1 



- rjf) sign n 



jJL — Z 



d/i. 



Then the general solution receives the following form 



X{z) \m{z) + Eoip{z) + E^z - ^!^z 



-WoA,T{z) + C,z. 



Prom this general solution we can find the function M{z) 

M{z) = -E^z - Eoip{r]o, z)ch — + — — ^ + -wqAi—- + — - 

770 2 3 X[z) X{z) 



(6.11) 



Let us remove a pole at the solution (6.11) at infinitely remote point. Let us notice, 
that the function 

at 2; — >■ 00 has the following asymptotic 



z ^ 00. 



Considering, that function T{z) has exactly the same asymptotic at 

Ci — {Eoo — )Aoo- 



(6.12) 



28 



8 Coefficients of discrete and continuous spectra 



Now we will remove poles at the solution (6.11) at the points ir^o- Let us allocate in 
the right part of the solution (6.11) members containing the polar singularity at the point 
z = rjo. In the point vicinity z = rjo taking into account equality X{r]o) = it is carried 
out the following expansion 



M{z) 



z + r)o 



+ 



Wo , {l/S)woAiT{z)+r]o{E^-woAi/2)X, 



770 A'(77o) + (l/2!)A"(77o)(-2-?7o) + --- 

Prom here it is visible, that for pole elimination at the point 2; = 770 it is necessary to 
equate to zero expression in a square bracket, calculated at 2; = rjQ. Then we receive, that 

EoX'{Vo){vI - Vl) ch(wo/?7o) - ?7o^ooAoo 



WqAi 



therefore 



2£; ch — = ^o^i[2^(^o)/3 - Aoor/o] + 2r)oX^E^ 
° Vo A'(77o)(77? -?7o) 



The continuous spectra coefficient E{r]) we find from formulas (6.11) and (6.12) 



27rz ch^(r72 
V 



ril) 



(7.1) 
(7.2) 

(7.3) 



Difference M~^{ri) — M {rj) from the formula (7.3) we will find with the help formulas 
of the general solution (6.12). As result we receive the following expression 



E{ri) 



1 



27ri ch — (772 — 77^) 



-WqAi 



X+{rj) X-{ri) 



WqA 



-)XooV 



X+{ri) X-{ri) 



(7.4) 



Let us substitute equalities (7.2) and (7.4) in the equation on the field (6.4). We 
receive the following equation 

2T(r/o)/3 - Aoo77o Ac 



A E 



where 



\X^ 



2770 



1] df] 



'-J1 + IJ2) = 1, (7.5) 
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J2 



T^irj) T {rj) \ dr] 



27ri y I A+(?7) A~(r/) j rf — rjl 



-1 



Integrals Ji and J2 from the equation (7.5) can be calculated analytically by means of 
the theory of residues and contour integration. For the first integral we have 



Ji 



2m 



1 



1 



rjdrj 

A+(?7) \~ {rj) \ rf — rjl 



Res + Res + Res + Res + Res 9x x / x ■ 

00 ?7i -?7i r)o -T?o J {z'^ — r]i)X\z) 



Let us notice, that 



Res 



Res 



^ {z-^-TjDXiz) Aoo' ±m {z-^-viDXiz) 2X{rj,y 

Res ^ = ^2 

Hence, the integral is equal 



1 1 
Ji = -— + ^ + 



Aoo ' Ai ' (^o-^i)A'M' 



where 



Ai = A(7^i) = 1 - -. 

2^0 



In the same way we calculate the second integral 



J2 



1 f T2{v) 



-J 



2TriJ rf -r]lVX+{r]) A" (77) J 
-1 



Res + Res + Res + Res + Res 

L 00 ?7i -771 ?jo -VO 



drj ■- 
T{z) 



Let us notice, that 



Res 



0, Res 



z'^ — r)j)X{z) ' 



Res 



T(z) 



±- 



±7o {z'-r,l)X{z) {vi-V^VoY 
Hence, the integral is equal 

_ rfa)-r(-7yi) T(r/o) - T(-r7o) 
■^2 — — ; r 



2r;iAi 



(^0 - vD^'ivo) 
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Let us notice, that T{r]i) — T{—r]i) — — . Thus, it is definitively received 

1 2r(r?o) 
2cX,^ (rji-ril)X'{rioy 

By means of these integrals the condition on the field (7.5) transforms to the form 

because Aoo = Ai + l/3c. 

Hence, from last equation we find the Drude amplitude 

Eoo = (7.6) 

Ann 



Now the condition (7.2) will be transformed to the form 

Ai?7o + woAi{lT{r]o) - IXooVo) 



^ _ ■•^■/u ' -~U"±V3- v/u/ 2—^ '"/ (7 7) 



9 Integral condition on distribution function 



Let us notice, that at transition through the positive part of the cut (0, 1) functions 
T{z) and X{z) make the jump, differing only by signs. Really, the formula for T{z) we 
will present in the kind 



or 



Tiz) 



z 

'2c 



fi — z IJ, + z 



1 1 
+ 



djji = - 

c 



2 2 "^Z^' 

/x^ — z"^ 



fjL — Z IjL + Z 



z ffi{fi^-r]l) 



//^ — z'^ 



This integral is easy for calculating in an explicit form. In a complex plane with a cut 
function T{z) is calculated by the formula 



f)ln(l-i) 

and on the cut this integral is calculated by formula 



zeC, 



-1 <r) < +1. 
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Now from the Sohkotsky formula for a difference of boundary values we receive, that 
at < 77 < 1 

X+(ri) = X{ri) ± ^Vivl - v'^), T+iv) = T{ri) ± -^viv'^ - vl) sign rj. 
Now it is easy to find, that 

T^v)X-{rj) - T-{r,)X^r)) = [T{r]) + X{r,) sign rj] ■ ^77(77^ - rjl). 



We enter the integral 



1 [ rf — r]l 



2c J r] — z 



drj. 



and we denote Ti(77) = T{r]) + X{ri) sign r]. It is clear that Ti(r/) is odd function, besides 



l + 2ijToi-r]), 7]>0, 
-{l-2riTo{ri)), 77 < 0. 



Substituting last two equalities in (7.3), we receive the formula for calculations of 
coefficient of continuous spectrum 



1 



2c ch (wo/rj) 



z^Ar 2r]T,{7])-3X 



+ 



Clearly that in complex plane the integral Tq{z) is calculated under the formula 



To{z) = 



2c 



l2 



+ z + {z^ -r]l 



and at < 77 < 1 it is calculated by formula 



2c 



^ + 77 + (r/' - ril) In - 1 



By means of this function we will present dispersion function in the form 

X{z) ^l- zn{z) + zTQ{-z), 
the function T{z) we express by this integral also 

T{z)^z%{z) + zT^{-z). 
The sum of two last expressions is equal 



X{z)+T{z) = l + 2zTo{-z). 
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We will notice, that the integral T{—z) is not singul on the cut < 77 < 1. The sum 
A(77) + T{r]) sign 77 on the cut of < 77 < 1 it is calculated in an exphcit form without 
quadratures 

X{ri) + T{ri) sign 77 = 1 + 2riTo{-ri) = 



r]-2r]^ + 277(77^ - rjl) In + ij 

Now we consider integral boundary condition (3.6). We rewrite it in the form 

1 

1 — CKp Ai 



CKp 36 



(8.2) 



Here, as it is easy to see. 



n[—a, ji) = ji 



F{rjo, /x)e"'°/^° + F(-77o, /x)e-"°/''°] + 



Wo Wo 
1 

+2 [ F{r],ii)E{r])e'^°/'^dr}. 



(8.3) 



Let us substitute (7.4) in expansion (7.5) for the function h{x,iJ,). We receive the 
following equation 

1 



E^ 



+ £^o7n(77o) + 2 / e'^°''^m{r])E{7])dr) = 

oD 



1 — CCp WqAi 

36~' 



Oir 



(8.4) 



where 



m(77o) = 6"'°/"° 7710 (770) + e-"'°/''°7no(-?7o). 
Besides, in the equation (8.4) following designations are accepted 
1 1 
mo{±rjo)=J (^/i^ -'^fi^F{±r]o,ii)dii, m{r)) ^ J (^fi^ - '^f?jF{r), ii)dii. 


Let us substitute in the equation (8.4) the Debayc amplitude £"0, defined by equation 
(7.7), and continuous spectrum coefficient £'(77), defined by equality (8.1). We receive the 
equation 



In the equation (8.5) are entered the following designations 

mo{vo)e'"°^'"' + m(-77o)e-"'o/'?o 



A'(?7o)(^o -^f)ch(wo/77o) 
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m{rj)e 



wo/r) 



X+{rj)\-{rj) ch(wo/?7)' 

-1 

Now from equation (8.5) we find the constant Ai 

^1 (set - ^0^0 + lo) 



rjTi{r))/3{r})dr). 



ZqAi 



(8.6) 



The solving of this initial boundary problem on it is finished. 

Let us calculate integrals mo(±?7o) and m{r]) in an explicit form. Integrals mo(r]o) and 
mo(— ?7o) are easily calculated 



mo{-rio) = ivf-Vo) 
We calculate the integral m{r]) 



(-^ + .o) + (.o^-^r,o)ln(l-i ^ 
-^-r;o) + (r/o^ + ^r,o)ln(l + i-) 



V 



di^-2c{r)--)X{r))e+{r)). 



ix-r] 



Here ^+(77) is the characteristical function of interval (0, 1), i.e 

1, 0<r/<l, 

0, -1<77<0. 

Prom here we receive, that at < v < 1 



1 2 1 2 

m{vi) = {nf - ril){- -7])- {nf - rii){iif - -r]) ln(l + -) + 2{ri^ - c){ri - -), 



and at —1 < 77 < 



rl 2 1 

m{ri) = {ri^ - ril) ---q- (77^ - -77) ln(l - -) 
Lb 3 77 J 



Last two formulas it is possible to unite in one 

m{v) = iv' - Vl) il-V-iv'- lv)U{v)] + 2(77^ - c)(77 - ^)^+(77), 

where 



Uiv) 



, 1 + 77 

In -, 0<77<1, 

, 1 — 77 

In -, -1 < 77 < 0. 
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10 Energy absorption in layer 

The energy of an electromagnetic wave absorbed in the slab of degenerate plasma bi 
well-known formula [35] 

a 

g^lRe j j{x)E*{x)dx. (9.1) 

—a 

Here "asterisk" means complex conjugation, and j{x) is the density of a current, 

1 

-1 

In view of one-dimensionality of a problem the equation for electric field has the 
form dE/dx — Attq. All quantities have dependence from time exp{—iujt), i.e. E — 
E{x) exp{—iujt), and etc. The continuity equation for system a charge - current will be 
rewrite so: dj{x)/dx — iuiq{x) — 0. Prom last two equations we have 

dE , djix) 
lu— = 47r- ^ ' 



dx dx 

Integrating it equality and considering equality to current zero on border, we have 

3{x)^'-^(e{x)-Eo), (9.2) 

where Eq is the amplitude of an external field on border (real quantity). 
Subsituting (9.2) into (9.1), we find 

a 

Q = Re / i{E{x) - Eo)E*{x) dx. (9.3) 



Stt 

—a 

We notice, that EE* = lE'p > is the real quantity, therefore. 

Re {i(E - Eo)E*} = -Eo Re (iE*) = 

= -Eo Rc {i{Re E - ilm E)) = -Eq Im E{x). 
Hence, according to (9.3) we have 

a 

coE f 
Q ~ — - — / Im E[x) dx, 



Stt 



—a 



or, by dimensionless electric field. 



Q^-'^^lm f e{xi)dxi. (9.4) 
Stt J 
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Here the expression received above for electric field is used, which we will present in 
the form 



-1 



A. 



V 



(9.5) 



Coefficients of continuous and discrete spectra of expansion (9.5) have been calculated in 
the explicit form. 

Let us substitute in the formula (9.4) expression for the electric field (9.5). After 
integration it is received, that 

V — Z VO; 

Air 

where Qq is the dimensionless part of absorption, defined by formula 

Qo = -lmQi, (9.5') 

where 

1 

Q, = ^ + 2Eo^sh^ + - [ rjEin) sh ^drj. (9.6) 
Aoo Wo rjo Wo J rj 

-1 

The formula (9.5) represents quantity of absorption of electric field energy in the metal 
layer. It is defined by frequency and amplitude of an external field u and Eq, width of 
the layer oq, and by parametres wq, zq, e, Q, refiecting properties metal. 

We consider integral from (9.6) 



J = rjEirj) sh — 
J V 



dr] 



2c 



1 1 



\\ — woA\l2 f rj th.{wo/r))dr) wqAi f rj ih.(wo/7])Ti(r))dr) 



+ 



/ 



A+(77)A-(77) 



Then taking into account equality 

1 c 



we receive 



where 



I ^1 

(n) X-(n)\ 



X+{rj)X-{rj) inri{rj^-rjf)lX+{r}) X-{r}) 



J-{X.-^)J. + ^Jo, 



rf sh(-u;o/''?) 
2m J {rf — rjl) ch.{wo/r]) 



x+iv) X-{7j) 



dr], 



-1 



rfTi{ri) th(wo/r7) _ 1 
A+(r;)A-(r;) " c 



ri'^Ti{r])i\i{wo/vi) 

x+{v)x-{v) 



dr). 



(9.7) 
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Let us notice, that integral Ji it is possible to calculate analytically, and integral Jq 
analytically to calculate it is impossible, since function continuation Ti{z) from the real 
axis in a complex plane it is carried out by non-analytical function \z\{l + 2zTo{—\z\))/z. 

Function zero ch{wo/z) will be necessary for us 

The integral Ji has been calculated earher in our work [26]. By means of contour 
integration and the theory of residues it is found, that 

Ji- 



Res + Res + Res + Res + Res + \^ Res 

k=—oo 

where 

z"^ th{wo/z) 



x{z){z'-viy 

Calculating all residues from the previous equality, we receive 

+00 4 



Jl = -T^ + 



wq 2?7gth(wo/?7o) rii wp 1 ^ 4 



Qi = ^ + 2£;o^ sh ^ + - X^^) J, + ^Jo. (9.9) 



By means of this expression the formula (9.6) can be presented in the explicit form 

) — sh h — 

In the formula (9.9) Debaye amphtude £'0 is calculated according to (7.7), quantity 
Ai is calculated according to (8.6), and integrals Jq and Ji are calculated under formulas 
(9.7) and (9.8) accordingly. 



11 Conclusions 

In the present paper the linearized problem of plasma oscillations in layer (particularly, 
thin films) in external longitudinal alternating electric field is solved analytically. Specular 
- accomodative boundary conditions of electron reflection from the plasma boundaries are 
considered. Coefficients of continuous and discrete spectra of the problem are found, and 
electron distribution function on the plasma boundary and electric field are expressed in 
explicit form. 

Separation of variables leads to characteristic system of the equations. Its solution in 
space of the generalized functions allows to find the eigen solutions of initial system of 
equations of Boltzmann — Vlasov — Maxwell, correspond to continuous spectrum. 
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Then the discrete spectrum of this problem consisting of zero of the dispersion equa- 
tions is investigated. Such zeros are three. One zero is infinite remote point of complex 
plane. It is correspond to the eigen solution "Drude mode". Others two points of dis- 
crete spectrum are differing with signs two zero of dispersion function. These zero are 
correspond to the eigen solution named " Debay mode" . 

It is found out, that on plane of parametres of a problem {Q, e), where Q = w/cjp, e = 
u/upj there is such domain (its exterior is domain D~), such, that if a point (0,£) G 
Debay mode is absent. 

Under eigen solution of initial system its general solution is obtained. By means 
of boundary conditions in an explicit form expressions for coefficients of discrete and 
continuous spectra are found. Then in explicit form (without quadratures) absorption 
quantity of energy of electric field in slab of degenerate plasmas is found. 
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